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ABSTRACT 

We describe the H l component of the Southern Galactic Plane Survey 
(SGPS). The SGPS is a large-scale project to image at arcminute resolution 
the H I spectral line and 21 cm continuum emission in parts of the plane of the 
Milky Way. The survey covers Galactic longitudes 253° < I < 358° and latitudes 
|6| < 1?5 (SGPS I), plus a hrst quadrant extension covering 5° < / < 20° and 
|6| < 1?5 (SGPS II). The survey combines data from the Australia Telescope 
Gompact Array and the Parkes Radio Telescope for sensitivity to angular scales 
ranging from 2 arcminutes to several degrees. The combined data cover 325 deg^ 
and have an rms sensitivity of 1.6 K. Here we describe the H l observations and 
data reduction in detail, and present examples from the hnal data products. The 
complete dataset is publicly available through the Australia Telescope National 
Facility’s H l Surveys archive. This dataset provides an unprecedented view of 
the neutral component of interstellar hydrogen in the inner Milky Way. 
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1. Introduction 

The physics of the interstellar medium (ISM) is diverse and complicated. The cycle of 
stellar formation and death takes the ISM through a wide range of physical conditions, a 
number of which can be traced by the 21-cm hyperhne emission line produced by neutral 
hydrogen (H l). Tracing both structure and dynamics, H l is observed in emission and 
absorption, each probing different thermal states of the gas. Recent observations by Heiles 
& Troland (2003) and Dickey et ah (2003) show that H l observed in absorption traces gas 
with temperatures as cold as 15 K, whereas H l emission observations generally trace gas with 
temperatures in the range of few thousands (Dwarakanath et ah 2002). H l densities vary over 
a similarly broad range, from ~ 10“^ cm“^ in the diffuse warm neutral medium to ~ 10^ cm“^ 
in regions of the cold neutral medium. The Milky Way provides an ideal laboratory in which 
to study the neutral ISM at sensitivities and spatial scales that are unattainable in external 
galaxies. Despite many years of observations of the ISM, our understanding of the medium 
is still simplistic, relying heavily on a three-phase model (e.g. McKee & Ostriker 1977). We 
have yet to adequately model the variety of temperatures, densities and structures known 
to exist in the ISM. Our lack of understanding is not only from the theoretical standpoint, 
but also from the lack of observational data probing H l physics over a wide range of spatial 
scales. It is clear that there is a great deal to learn about the structure, dynamics and 
thermal state of the neutral ISM from high resolution surveys of Galactic H l. 

To reach this goal and develop a better understanding of the Galactic ISM the Interna¬ 
tional Galactic Plane Survey (IGPS) is systematically imaging H l in the disk of the Milky 
Way. The IGPS consists of three individual surveys, the Ganadian Galactic Plane Survey 
(GGPS; Taylor et ah 2003), the VLA Galactic Plane Survey (VGPS; Taylor et ah 2002), and 
the Southern Galactic Plane Survey (SGPS, described here). These surveys excel because of 
their exceptional spatial dynamic range, combining single dish and interferometer data for 
arcminute resolution imaging over areas of many degrees. 

The Southern Galactic Plane Survey (McGlure-Grifhths et ah 2001b) is a survey of the 
21-cm continuum and H I spectral line emission in the fourth and parts of the first and third 
quadrants of the disk of the Milky Way. The survey provides a complete H l and 21 cm full 
polarization (Stokes I, Q, U, & V) dataset of the area 253° < I < 358°, b = ±1?5 (SGPS I) 
and an extension into the hrst quadrant covering 5° < / < 20°, b = ±1?5 (SGPS II). The 
SGPS is comprised of two surveys: a low resolution (FWHM = 15') survey completed with 
the Parkes 64m radio telescope and a high resolution (FWHM = 2') interferometric survey 
completed with the Australia Telescope Gompact Array (ATGA). The data are combined to 
provide a dataset that is sensitive to more than two orders of magnitude in spatial scales. The 
combined data have an angular resolution of ~ 2 arcmin, spectral resolution of 0.8 km s“^ 
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and a sensitivity limit of ~ 1.6 K in the line. The main SGPS parameters are given in 
Table 7. 

The scientihc objectives of the SGPS are to provide a high-resolntion 21cm polarized 
continnnm and H l atlas of the Galactic plane with which to stndy the distribntion and 
dynamics of both the warm and cool components of the ISM over size scales ranging from 
parsecs to kiloparsecs (McGlnre-Grifhths et ah 2001b; McGlnre-Grifhths 2001). Scientihc 
results from the survey and its “test region” have been presented in a number of papers. An 
overview of the H l emission and absorption in the test region was published in McGlnre- 
Grifhths et ah (2001b). Other SGPS papers focusing on the test region data include an 
analysis of the H l emission spatial power spectrum (Dickey et al. 2001); a discussion of the 
abundance and temperature of cool H l clouds in the inner Galaxy (Dickey et al. 2003); an 
examination of extended polarized emission structure and rotation measures of point sources 
(Gaensler et al. 2001); and an analysis of rotation measure huctuations indicating an outer 
scale for turbulence of ~ 2 pc (Haverkorn et al. 2004). Using the SGPS as a whole we 
have discovered a new, outer spiral arm in the fourth quadrant of the Milky Way (McGlnre- 
Grifhths et al. 2004); probed the physics of individual H l supershells (McGlnre-Grifhths et al. 
2000, 2001a, 2003); discussed the Galactic distribution of H l supershells from a catalog of 
new supershells discovered in the SGPS (McGlure-Grifhths et al. 2002); and explored the 
properties of H l self-absorption clouds discovered in the SGPS (Kavars et al. 2003, 2004). 

Here we describe the SGPS in detail, focusing on the observing strategy and data re¬ 
duction techniques. In particular, we describe the H l data; the continuum and polarization 
data require vastly diherent reduction techniques and will be described in a separate pa¬ 
per (M. Haverkorn et al. 2004, in preparation). This paper is organized as follows: in §2 we 
describe the low-resolution Parkes multibeam survey, in particular the data reduction strate¬ 
gies developed explicitly for this project. In §3 we describe the ATGA observations and data 
reduction, especially where the reduction deviates from standard techniques. The combi¬ 
nation of our single dish and interferometer data is described in §4. Finally, the complete 
data products are described in §5 and representative images from the combined data cubes 
are presented in §5.2. The SGPS data are publicly available from the Australia Telescope 
National Facility (ATNF) H l Surveys archive.^ 


^http://www.atnf.csiro.au/research/HI/common/ 
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2. Parkes Survey 

The low resolution portion of the SGPS was observed with the Parkes Radio Telescope. 
This survey covers the area 253° < I < 20°, b = ±10° with an angular resolution of 15'. The 
Parkes survey was primarily designed to provide information about large scale structure not 
sampled in the ATCA survey. Expanding the latitude coverage of the Parkes survey to trace 
large-scale structures away from the Galactic plane provides additional information about 
the gas in a thicker disk and has led to the discovery of many new large H l shell structures 
(McGlure-Grifhths et ah 2002). Therefore, the Parkes survey is released as a separate data 
product. Here we describe the observational and data reduction techniques used for this 
survey. 


2.1. Observations 

Single dish observations were made with the 21cm multibeam system at prime focus on 
the Parkes 64m Radio Telescope near Parkes, NSW, Australia. This is a focal plane array 
with thirteen independent feeds, each with dual, cryogenically cooled, orthogonal linear 
polarization receivers with stable, low system temperatures {Tgys ~ 25 K). The feeds are 
packed in a hexagonal pattern with a separation of 29H between the centers of adjacent 
feeds (Staveley-Smith et al. 1996). At A = 21 cm the FWHM of the Parkes beam is 14f4, so 
the multibeam array does not fully sample the focal plane. For this project only the inner 
seven beams were used. 

Data were obtained by mapping “on-the-fly”, driving the telescope at a constant Galac¬ 
tic longitude through 3 degrees of Galactic latitude at a rate of 1.1 deg min“^. Data samples 
were written to disk every 5s. In order to maximize sky coverage and reduce redundant sam¬ 
ples the multibeam receiver platform was rotated to maintain an angle of 19H between the 
face of the hexagon and the scan direction. The platform rotation resnlts in seven parallel 
beam tracks, each separated by 9^5. To ensnre Nyquist sampling, interleaved scans were 
made by offsetting the receiver package by 2.5 beam spacings and scanning in the opposite 
direction. The offset scans reduced the spacing between adjacent tracks to 4T and ensured 
that independent feeds are responsible for adjacent tracks. Along the scan direction the data 
samples are spaced by ~ 5'5, which smears the beam slightly in the scan direction. In order 
to rednce the effects of variations in the system temperatnre with elevation, we made an 
effort to scan the telescope at a nearly constant zenith angle of approximately 30°. 

Observations were condncted dnring Eve sessions between December 1998 and November 
2000. The first set of observations on 15 and 16 December 1998 covered the ATGA survey 
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region of 253° < I < 358°, |6| < 1?5. We conducted repeat observations of this same region 
at three month intervals for a year in order to reduce the annual variations in H l spectra due 
to stray radiation (Kalberla et ah 1980) and to improve the sensitivity of the observations. 
Ultimately the ATCA survey region was observed 10 times. During subsequent observing 
runs in June 1999, September 1999, and March 2000 the latitude range was extended to 
cover |6| < 10°. Finally, between 28 October and 5 November 2000 the survey was extended 
through the Galactic Center to I = 20°. 

The data were recorded in frequency switched mode, switching in band between a cen¬ 
ter frequency of 1419.0 and 1422.125 MHz with every 5s cycle. We used the narrow-band 
multibeam correlator (Haynes et al. 1999), with a maximum bandwidth for each frequency 
of 8 MHz, with 2048 channels, giving a channel width of 3.9 kHz (An = 0.82 km s“^). After 
bandpass correction using frequency switched samples the usable bandwidth was approxi¬ 
mately 4.5 MHz. 


2.2. Bandpass and Brightness Temperature Calibration 

The lAU standard narrow line calibration regions S6 and S9 were observed daily for 
absolute brightness temperature and bandpass calibration. For these observations each beam 
of the array was pointed on source for a 60s integration. These regions were also observed 
near a zenith angle of 30° so as not to compromise calibration with a signihcantly different 
gain and system temperature from the survey helds. 

The bandpass shape is relatively stable over periods of days. However, there are vari¬ 
ations in the front-end gain versus frequency shape on the timescale of seconds to minutes 
that can largely be removed by frequency switching. To remove the bandpass ripples, each 
spectrum was divided by the previous frequency switched sample to create a quotient spec¬ 
trum. After this hrst reduction step, the shape of the quotient spectrum is stable to the 1-3% 
level over the course of a day. Next, the residual bandpass shape was corrected by dividing 
the quotient spectrum by a bandpass template. The bandpass template was determined 
daily for each beam and polarization from a time-averaged quotient spectrum on one of the 
standard line regions, S6 or S9. These spectra were normalized to a mean of one and then 
hand-£t by a series of Fourier components to create a bandpass template, /(zz). In general, 
the narrower line width of S6 (An = 10.6 km s“^ for S6, versus An = 15.8 km s“^ for S9) 
resulted in a better £t for the template. 

In standard processing of H l single dish spectra the antenna temperature is calculated 
as Ta = T'y^{Tobs - Tref)T~^), where is the system temperature determined online during 
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the integration, Tabs is the observed spectrum and T^e/ is the subsequent frequency-switched 
spectrum. The application of this technique, with some polynomial bandpass htting, is 
relatively successful at removing ripples from the bandpass shape, but it also removes any 
continuum emission. For the SGPS we are interested in the continuum emission as well as 
the H I line so we employed a non-standard technique. For each spectrum we calculated 
the mean value, Tbas, of the reference spectrum over the area of the line, which should be 
the sum of the sky continuum temperature and the actual system temperature, Tgyg. For 
Tsys we used an average system temperature calculated on the standard line regions, which 
we apply to an entire day’s observations. This is a generally valid assumption because the 
system temperature of the multibeam is extremely stable. The bandpass corrected spectra 
were then calculated as: 


Ta{iy) = T, 


Tobsi^) 


/(O 

where is a smoothed version of the frequency switched, or reference, spectrum. This 

technique differs from the standard single dish H l bandpass calibration because it only 
subtracts the absolute system temperature, not the sum of the system temperature and 
the continuum emission. The resultant spectra contain the observed continuum emission 
and have baselines that are typically flat to less than IK. They are limited, however, by the 
accuracy in our determination of Tgys- The value of Tgjs contains not only the receiving system 
temperature, but also the ground radiation temperature, which is elevation dependent and 
poorly known at 21 cm the for Parkes telescope. Because of the difficulty in estimating total 
power continuum levels with a single dish, Tgys is only accurate to ~ 10 — 15%. Fortunately, 
for data products from which the continuum emission is subtracted later in the analysis this 
problem is alleviated. 


After bandpass correction, the integral over part of the line on the standard regions was 
used to determine a calibration scale factor, G, to convert antenna temperature to absolute 
brightness temperature, Tb = C Ta. The calibration factor was calculated for each beam and 
each polarization by comparing the measured line integral on the standard line regions to 
the published Williams (1973) values for the line integral, such that C = hrue/Iobs- For S6 
the line integral was calculated over the velocity range —5.8 km s“^ < n < -1-4.8 km s“^ 
and is Itme = 299 K km s“^. The S9 line integral over the velocity range -1-1.05 km s“^ < 
V < -1-14.75 km s“^ is Itme = 953 K km s“^. Using the line integral to scale the measured 
antenna temperature we hnd that the peak brightness temperature for S6 is T^ax = 57±2 K 
and Tmax = 90 ± 3 K for S9. For comparison, another recent Parkes multibeam survey with 
a similar observational setup found Tmax = 83 K for S9 (Briins et ah 2003). The discrepancy 
may be due to the difference between the velocity range used for the line integral in Williams 
(1973), from which our calibration was derived and Kalberla et ah (1982), from which the 
Briins et ah (2003) calibration was derived. 
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Finally, the spectra are shifted to the Local Standard of Rest (LSR) reference frame. 
The correction is calculated using standard techniques and is applied as a phase shift in the 
Fourier domain. 


2.3. Imaging 

The fully calibrated and velocity shifted Parkes spectra were imaged using Gridzilla, a 
gridding tool created for use with Parkes multibeam data and found in the ATNF subset of 
the aips++ (now casa) package. The gridding algorithm is described in detail by Barnes et 
ah (2001). The SGPS data were gridded using a weighted median technique with a cellsize 
of 4', a Gaussian smoothing kernel of 16', and a smoothing kernel cutoff radius of 10'. The 
Gaussian smoothing kernel and cutoff radius reduce the effective angular resolution of the 
gridded data to ~ 16', which is slightly larger than the lA'A FWHM of the Parkes beam at 
1420 MHz. The effective resolution is estimated from Gaussian fits to point sources in the 
Parkes continuum images and to a beam map made on the calibration source PKS B1934- 
638. The imaged data were imported into the MIRIAD data reduction package (Sault & 
Killeen 2004) for further analysis and continuum subtraction. Image domain continuum 
subtraction was performed using a first order polynomial £t to the off-line channels. 


3. Australia Telescope Compact Array Survey 

The high resolution portion of the SGPS was completed with the Australia Telescope 
Gompact Array (ATGA) synthesis imaging telescope. The ATGA survey is divided into two 
parts: SGPS I, covering the area 253° < / < 358°, |&| < 1?5 and SGPS II, covering the area 
5° < / < 20°, |6| < 1?5. The Galactic center area from 358° < / < 5° has also been observed, 
but will be presented separately because it has a different observing and reduction strategy. 
Here we describe the observational and data reduction strategies used for both SGPS I and 
H. 


3.1. Observations 

The observations for SGPS I were conducted between December 1998 and August 2000, 
with re-observation of several areas in June 2001. SGPS H was observed between June and 
October 2002. The observational strategies are slightly different for the two data sets. 

The ATGA is a six element interferometer near Narrabri, New South Wales, Australia. 



It consists of five 22 m antennas that can be positioned along on a 3 km east-west track 
or a perpendicular 200 m north-south track and a sixth antenna hxed a further 3 km west 
of the end of the east-west track, for a maximum baseline of 6 km. The antennas can be 
arranged in a number of conhgurations. The SGPS takes advantage of the ATCA’s most 
compact conhgurations with maximum baselines (excluding the 6 km antenna) between 210 
m and 750 m to maximize sensitivity to large scale structures. The array conhgurations 
used for SGPS I were: 210, 375, 750A, 750B, 750G, and 750D. The array conhgurations were 
chosen to ensure optimum sampling of the inner u-v plane and minimize redundant baselines. 
All baselines, in intervals of 15.5 m, from 31 m to 291 m, except 199m, are sampled. The 
long baselines of the 750-m arrays are poorly sampled and limit our practical resolution 
to ~ 2'. Before phase II of the SGPS the ATGA was upgraded to include a short north- 
south track and to add new stations on the east-west track that improve the coverage of the 
compact conhgurations. Because the SGPS II extends as far north as 5 = —11°, the new 
north-south baselines were utilized to improve the u-v coverage and to reduce the ehects of 
shadowing at low-elevations. The SGPS II was observed with the new EW352 and EW367 
conhgurations, which give nearly complete coverage to baselines of 352 m in intervals of 
15.5 m, and the north-south/east-west hybrid arrays, 1175 and 11168. Even with the hybrid 
arrays, the synthesized beam for the more northerly of the SGPS II helds is quite elliptical 
(200" X 130"). 

In order to cover an area much larger than the primary beam of the ATGA antennas 
(Gpb = \/D = 33'), the survey was conducted as a mosaic of many pointings. The pointings 
were arranged in a hexagonal pattern with Nyquist spacing of 19 arcmin between adjacent 
pointings. The theoretical rms sensitivity for Nyquist sampled pointings is uniform to within 
< 1% across the entire observed area. In the SGPS I there are a total of 2212 pointings 
comprising twenty helds of 105 pointings and one held of 112 pointings. In the SGPS II 
there are 315 pointings divided in three helds of 105 pointings each. The helds composed 
of 105 pointings cover 5?5 x 3°, centered on a Galactic latitude of 6 = 0°. Each pointing 
was observed approximately 40 times with 30 second integrations for a total integration time 
of at least 20 minutes. The individual 30 second integrations are distributed evenly with 
hour angle. The exact u-v coverage varies from pointing to pointing, but is mostly constant 
within any given held. An example of the u-v coverage for one pointing near the center of 
the held at I = 295° is shown in Figure 1. The u-v coverage for this pointing is typical of 
the SGPS I helds. Similarly, an example of the u-v coverage for a pointing near the center 
of the held at I = 12° is shown in Figure 2. This coverage is typical of the SGPS II helds, 
although the coverage is much more circular for pointings near I = 7° and more elliptical for 
pointings near I = 17°. The u-v coverage shown in Figure 2 dihers from the coverage shown 
in Figure 1 because of the use of the north-south hybrid arrays in SGPS II. 
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The ATCA has linear feeds, receiving two orthogonal linear polarizations, X and Y, 
and is capable of observing two frequencies simultaneously. All data were recorded in a 
spectral line mode with 1024 channels across a 4 MHz bandwidth centered at 1420 MHz, and 
simultaneously a continuum mode with 32 channels across a 128 MHz bandwidth centered 
at 1384 MHz. In continuum mode the full polarization products, XX, YY, XY, and YX 
are recorded, whereas only the autocorrelations XX and YY are recorded in spectral line 
mode. Only the narrowband data are described here. 


3.2. Calibration and Imaging 

Data editing, calibration and imaging were conducted in the MIRIAD data reduction 
package using standard techniques (Sault & Killeen 2004). The primary flux calibrator for 
the southern hemisphere, PKS B1934-638, was observed at least once per day and was used 
for bandpass and absolute flux calibration of all data, assuming a flux at 1420 MHz of 
14.86 Jy (Reynolds 1994). Observations of the secondary calibration sources PKS B0823- 
500, PKS B1105-680, PKS B1236-684, and PKS B1833-210 were used to solve for the time 
varying gain and delay corrections. One secondary calibrator was observed approximately 
every hour during the course of observations. The source PKS B0823-500 was observed for 
all helds between I = 255° and I = 290°. PKS 1105-680 was observed for all helds from 
I = 290° to I = 300°, and PKS 1236-684 was observed for all helds between I = 300° and 
/ = 320°. For helds between I = 320° and I = 355° PKS 1934-638 was used as both a primary 
and a secondary calibrator. Finally, PKS B1833-210 was used as the secondary calibrator 
for the three helds in SGPS H. In all cases, the gain solutions from the nearest secondary 
were copied to the SGPS held data. For helds hanked by two secondary calibrators the two 
solutions were merged when copied to the SGPS held. 

The individual pointings of a held were linearly combined and imaged using a standard 
grid-and-FFT scheme. The mosaicing process uses the joint approach where dirty images of 
each pointing are linearly combined and then jointly deconvolved (Sault et ah 1996). This 
approach can oher a better deconvolution because information away from the pointing center 
is included in the deconvolution process. Additionally, this method takes best advantage of 
the mosaicing technique’s ability to recover information on angular scales larger than the 
scale determined by the minimum baseline separation, \/dmm ~ 23', where dmin is the 
minimum baseline. Jointly imaging and deconvolving the pointings increases the maximum 
angular scale imaged to A/((Jmin — D/2) ~ 36', where D is the diameter of a single antenna. 

Super-uniform weighting was used for all H l cubes. Like uniform weighting, super¬ 
uniform weighting minimizes sidelobe levels to improve the dynamic range and sensitivity to 
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extended structures. However, uniform weighting reverts to natural weighting if the held of 
view is much larger than the primary beam, as is the case for large mosaics. Super-uniform 
weighting decouples the weighting from the held size and is therefore more successful than 
uniform weighting on large mosaics (Sault & Killeen 2004). 

Two cubes of each imaged region were made: one containing continuum emission and 
one for which the continuum had been subtracted. For the continuum subtracted H l cubes, 
the continuum emission was subtracted from the calibrated u-v data prior to imaging (Sault 
1994). These data were imaged without the 6 km antenna because of the poor sensitivity 
of the long baselines. For the continuum subtracted cubes, the pointings from two adjacent 
helds were jointly imaged to produce cubes covering 11° x 3° {I x b). For the H l cubes 
containing continuum emission the 6 km baselines were retained for use with H l absorption 
experiments. All cubes were then deconvolved using a maximum entropy deconvolution 
algorithm (Sault et ah 1996). The continuum-subtracted cubes for SGPS I were restored 
with a circular synthesized beam of 130". Because of their more northernly declinations, 
the SGPS II cubes are restored with slightly elliptical beams, 160" x 110" for field gOlO and 
200" X 130" for held g015. The SGPS I and II cubes with continuum emission and the 6 km 
antenna were restored with varying beam sizes, given in Table 4. 


4. Single Dish and Interferometer Data Combination 

Although mosaicing recovers information on larger angular scales than normal interfer¬ 
ometric observations by reducing the effective shortest baseline, the ATGA images alone are 
not sensitive to angular scales larger than ~ 36 arcmin. In order to recover information on 
scales as large as the image size, we have combined the ATGA data with the Parkes data of 
the same region, which continuously sample the u-v plane between zero baseline separation 
and the 64 m diameter of the dish. There are a number of methods for combining single dish 
and interferometric data. The merits of the various methods have been discussed extensively 
by others, (e.g. Stanimirovic 2002), so we will not repeat the discussion. For the SGPS we 
combine the data in the Fourier domain after deconvolution, as implemented in the MIRIAD 
task IMMERGE. We chose this combination method for several reasons, specihcally that it 
less sensitive than other methods to errors in the model for the Parkes beam, is not overly 
computationally intensive and produces reliable results. The ATGA and Parkes data are 
deconvolved separately, Fourier transformed, reweighted, linearly combined and then inverse 
Fourier transformed. If Vint{k) is the Fourier transform of the deconvolved ATGA mosaic and 
Vsd{k) is the Fourier transform of the deconvolved Parkes image, then the Fourier transform 
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of the combined image is given by 


Vcomb{k) = io'{k)Vint{k) + fcal^''{k) Vsd{k). (2) 

The weighting functions, uj'{k) and uj"{k), are dehned such that u'{k) +uj"{k) is a Gaussian 
whose full-width half-max is the same as the synthesized ATCA beam. In this way the 
ATCA data are downweighted on short baselines and the Parkes data are downweighted 
on the longer baselines. The calibration factor, fcai, scales the Parkes data to match the 
brightness temperature scale of the ATCA data. This factor is determined by comparing the 
ATCA and Parkes datasets at every pixel and frequency in the range of overlapping spatial 
frequencies, which for an ATCA mosaic and a Parkes image is ~ 120A to ~ 190A, where 
A = 21 cm. For the SGPS we found that fcai was equivalent to one to within 2%, which 
adds conhdence to the calibration of the individual Parkes and ATCA datasets. The hnal 
combined H l cubes are regridded to Galactic coordinates with a Cartesian projection. 


5. Resultant Data Products 

There are three resultant H l data products: the Parkes H l cubes, ATCA-|-Parkes 
continuum subtracted H l cubes and ATCA-|-Parkes H l cubes with continuum emission. All 
three data products are available for viewing and downloading at the ATNF H l Surveys 
Archive at http://www.atnf.csiro.au/research/HI/common. 


5.1. Parkes Data Cubes 

Because the Parkes survey covers a much larger area than the ATCA survey these data 
are made available separately. The Parkes data are presented as twelve cubes, most covering 
17° X 20° (/ X b). Adjacent cubes overlap by about 4 degrees in longitude to allow for 
convenient subsequent combination. The center position, size and sensitivity of the Parkes 
cubes are given in Table 2. The hnal cubes have an angular resolution of 15' and a velocity 
resolution of 0.82 km s“^. The rms noise in line-free channels is approximately 180 mK for 
the regions 1?5 < |&| < 10°, and approximately 60 mK in the region \b\ < 1?5, which was 
observed many times. However, stray radiation from the far sidelobes can cause spurious 
emission with amplitudes as high as ~ IK, so low brightness temperature regions should be 
viewed with caution. As an example of the data quality we display the entire Parkes dataset 
at a velocity of u = 34.4 km s“^ in Figures 3. These images show the variety of size scales 
represented in the Parkes data alone. 



12 


5.2. Combined Parkes and ATCA Cubes 

A description and samples of the final combined Parkes and ATCA cubes are presented 
here. The data are available from the ATNF H l Surveys Archive. The cube sizes, cen¬ 
ter positions, line-free rms noise per channel and synthesized beam sizes for the continuum 
subtracted data and data containing continuum emission are given in Tables 3 and 4, re¬ 
spectively. All cubes have 40 arcsec pixels and a channel separation, 5v = 0.82 km s“^. The 
velocity range of the imaged cubes differs from cube to cube. These ranges were chosen to 
include the bulk of the Galactic emission and minimize cube size for download. The mean 
rms in the line-free channels is 1.6 K for both the SGPS I and II continuum subtracted 
cubes. The mean rms for the cubes containing continuum emission is 1.8 K. The positions 
and dimensions of the combined continuum subtracted cubes are shown in Figure 4. 

Because of the large volume of data in the SGPS it is impractical to display all channels 
of the data cubes here. Instead we have chosen to present only longitude-velocity {l-v) 
diagrams at 6 = 0° and an example velocity channel from one of the cubes. The l-v diagrams 
for SGPS I and II are shown in Figures 5 and 6, respectively. There are approximately 
2900 independent spectra in the SGPS I l-v diagram at 6 = 0° and 300 in the SGPS II 
l-v diagram at 6 = 0°. Figure 7 is a single velocity channel at v = 75.05 km s“^ from the 
continuum subtracted Parkes-|-ATGA g268 cube, compared with the same velocity channel 
in the Parkes data alone. This figure shows the dramatic increase in resolution of the SGPS 
combined image compared with just the single dish data. Figure 7 also demonstrates that 
the combined data accurately represent the large scale structure present in the single dish 
data. 


5.3. Data Quality and Artifacts 

We have taken care to reduce the number of artifacts present in the SGPS data, how¬ 
ever there are some artifacts that are unavoidable. These are, in general, artifacts related to 
gain calibration of the Parkes multibeam and limited u-v coverage in the ATGA data. The 
Parkes data have some residual longitudinal striping caused by inter-beam gain variations 
in the multibeam and errors in the estimation of the system temperature, Tgys, as explained 
in §2.2. For the continuum subtracted cubes, errors due to Tgys are subtracted, leaving only 
gain variations. These variations are least evident in the region \b\ < 1?5, which was observed 
many times, but are apparent in the higher latitude data. For the cubes containing contin¬ 
uum emission, errors in the estimation of Tsyg effect the reliability of brightness temperature 
values on the largest angular scales. Once again, this problem is mitigated for data in the 
region |6| < 1?5 where ten independently calibrated datasets were averaged together. 
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The Parkes cubes also contain some artifacts related to saturation. There are four 
positions along the Galactic plane where strong continuum emission saturated the data 
samplers at Parkes. At these positions the observed fluxes are unreliable so we have blanked 
the images. These appear as white spots in Figure 3. 

The ATCA data suffer from artifacts resulting from a limited number of snapshots and 
incomplete u-v coverage. Because the interferometric data contain a small number (usually 
~ 40) of short duration hour-angle cuts these data have a dynamic range of about 150:1. 
For Galactic plane H l emission the measured brightness temperature is generally less than 
150 K so the images are sensitivity limited, not dynamic range limited. However, towards 
the numerous bright continuum sources in the Galactic plane, the dynamic range is the 
fundamental limitation. These continuum sources appear in absorption in the continuum 
subtracted data where they can limit the images as much as in the non-continuum subtracted 
data. The limited u-v coverage also leads to a number of radial and azimuthal sidelobe 
artifacts around bright continuum sources that are not fully removed by the deconvolution 
process. The SGPS I data suffer from effects related to the fact that all baselines are multiples 
of 15.3 m. The result is a grating lobe of radius 0?8 around strong continuum sources. The 
addition of the Parkes data to the interferometric images reduces the effect of the grating 
lobe but it remains at a level of about 5% of the continuum peak. This artifact is further 
reduced in the SGPS II cubes, which contain north-south hybrid arrays with a wider variety 
of baseline spacings. 


6. Discussion 

SGPS data can be applied to a number of scientihc applications that demand high- 
resolution H I imaging of large areas. In the Galactic ISM many structures are large and have 
complex interrelationships, which are not apparent from more piecemeal observations. It is in 
these studies that the SGPS can excel, revealing hne detail that would not have been observed 
in the large single dish surveys previously available. However, a realistic 3-dimensional 
picture of the ISM is only possible if the large-scale flux density is measured, which is not 
measured in interferometric surveys. The large-scale flux is essential for producing a true 
estimate of the H l mass and dynamics of the Galaxy and understanding the nature of 
ISM structures. For example, there are some features which are principally perturbations 
in velocity space caused by turbulence whereas there are others that have massive density 
variations (e.g. shells with voids). Purely interferometric surveys, while showing details well, 
can fail to distinguish between these two types of features. 

This paper gives a description of the three data products of the SGPS surveys of the 
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H I line: Parkes continuum subtracted H I cubes with a resolution of 15', combined Parkes 
and ATCA continuum subtracted H l cubes with a resolution of 2'.2 and combined Parkes 
and ATCA cubes with a resolution of ~ 1'6 containing continuum emission. Each of these 
data products was produced for use in different general scientihc purposes. The Parkes H l 
data, though low resolution, cover an extended latitude range. These data are most useful 
for large-scale panoramic views of the H l, tracing Galactic structure and studies of discrete 
structures that extend beyond the thin disk of the Galactic plane. The combined Parkes 
and ATCA continuum subtracted cubes provide high-resolution H l images with maximum 
surface brightness sensitivity. These data are best-suited for high-resolution examinations of 
H I emission structures. The combination of the ATCA data with Parkes data ensures that 
these cubes can be used to explore the ISM on multi-degree scales. Finally, the combined 
Parkes ATCA data containing continuum emission provide accurately calibrated data at 
a slightly higher angular resolution. These data were produced specihcally for use in H l 
absorption studies, which require good knowledge of the continuum flux as well as the H l 
flux in order to accurately measure optical depth. 


7. Conclusions 

The Southern Galactic Plane Survey (SGPS) is a major new data resource for Galactic 
H I astronomy, providing an unprecedented view of H l in the inner Galaxy in the southern 
hemisphere. As part of the International Galactic Plane Survey, the SGPS contributes to 
nearly continuous coverage of the 1st, 2nd, and 4th quadrants of the plane of the Galaxy. 
We have presented the H l data products which are now publicly available via the ATNF H l 
Surveys archive. The H l data products include Parkes H l cubes with 15 arcmin resolution 
and ~ 100 mK rms sensitivity covering the area 253° < / < 20° and |&| < 10°, continuum 
subtracted combined ATGA and Parkes H l cubes with 2 arcmin resolution and ~ 1.6 K 
rms sensitivity covering the regions 253° < / < 358°, |6| < 1?5 (SGPS I) and 5° < / < 20°, 
\b\ < 1?5 (SGPS II), and combined ATGA and Parkes H l cubes of the same regions including 
the continuum emission. 


The ATGA and the Parkes Radio Telescope are part of the Australia Telescope which 
is funded by the Gommonwealth of Australia for operation as a National Facility managed 
by GSIRO. This research was supported by NSF grants AST-9732695 and AST-0307603 to 
the University of Minnesota and AST-0307358 to Harvard University. Additional support 
was provided through a NASA Graduate Student Researchers Program (GSRP) Fellowship 
to N. M. M.-G. while at the University of Minnesota. We would like to thank the staff 
of the Australia Telescope National Facility for their support of this project, especially M. 
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Table 1. Basic parameters for the SGPS H l survey 


Dataset 

Angular Resolution 

Spectral Resolution 

Longitude coverage 

Latitude coverage 

1-cr Sensitivity 

Parkes SGPS 

15' 

0.8 km s“^ 

253° < ; < 20° 

|5| < 10° 

180 mK (60 mK'^) 

ATCA + Parkes SGPS I 

212 

0.8 km s“^ 

253° < 1 < 358° 

\b\ < 1?5 

~ 1.6 K 

ATGA + Parkes SGPS II 

~ 313 

0.8 km s“^ 

5° <l < 20° 

6 < 1?5 

~ 1.6 K 


The rms in the latitude range |5| < 1?5. 



Table 2. Cube properties for the Parkes H l Cubes. Column headings are as follows: Cube name; center Galactic 
longitude, /; center Galactic latitude, 6; longitude pixel size, 61, latitude pixel size, 66] number of pixels in /, npix;; 
number of pixels in b, npix^; longitude coverage. A/; latitude coverage, A6; minimum velocity, Vmin] maximum 

velocity, Vmax] velocity channel width, 6v; and cube rms. 


cube 

1 

(°) 

b 

n 

51 

C) 

npix; 

5b 

(') 

npixfc 

Al 

n 

Ab 

n 

'^min 

(km s“^) 

^max 

(km s“^) 

5v 

(km s“^) 

rms 

(mK) 

p268 

264.5 

0.0 

4 

345 

4 

305 

23.00 

20.33 

-149.8 

220.4 

0.82 

130 

p278 

277.5 

0.0 

4 

255 

4 

305 

17.00 

20.33 

-149.8 

250.0 

0.82 

130 

p288 

287.8 

0.0 

4 

255 

4 

305 

17.00 

20.33 

-180.3 

250.0 

0.82 

140 

p298 

297.5 

0.0 

4 

220 

4 

300 

14.67 

20.0 

-200.1 

199.7 

0.82 

200 

p308 

307.5 

0.0 

4 

220 

4 

300 

14.67 

20.0 

-200.1 

199.7 

0.82 

210 

p318 

317.5 

0.0 

4 

225 

4 

300 

15.00 

20.0 

-350.2 

349.8 

0.82 

200 

p328 

327.5 

0.0 

4 

225 

4 

300 

15.00 

20.0 

-200.1 

250.0 

0.82 

200 

p338 

337.5 

0.0 

4 

220 

4 

300 

14.67 

20.0 

-299.9 

300.3 

0.82 

230 

p348 

347.5 

0.0 

4 

225 

4 

300 

15.00 

20.0 

-299.9 

199.7 

0.82 

190 

p358 

357.5 

0.0 

4 

225 

4 

300 

15.00 

20.0 

-299.9 

199.7 

0.82 

180 

p003 

3.0 

0.0 

4 

249 

4 

300 

16.60 

20.0 

-399.6 

381.1 

0.82 

170 

p015 

15.0 

0.0 

4 

225 

4 

300 

15.00 

20.0 

-149.8 

300.3 

0.82 

220 



Table 3. Properties of the ATCA and Parkes combined H l continuum subtracted cubes. Column headings are as 
follows: Cube name; center Galactic longitude, /; center Galactic latitude, b; longitude pixel size, SI, latitude pixel 
size, Sb] number of pixels in /, npix;; number of pixels in b, npix^; longitude coverage, Al; latitude coverage, Ab; 
minimum velocity, Vmin] maximum velocity, Vmax', velocity channel width, Sv, synthesized beam size, 9; and the cube 

per channel rms. 


cube 

1 

n 

b 

n 

51 

(") 

npix/ 

5h 

n 

npixfe 

Al 

n 

Ab 

n 

'^min 

(km s“^) 

'^max 

(km s“^) 

Sv 

(km s“^) 

9 

(" X ") 

rms 

(K) 

g258 

257.9 

0.0 

40 

990 

40 

260 

11.00 

2.89 

-112.1 

216.8 

0.82 

130 X 130 

1.4 

g268 

267.9 

0.0 

40 

990 

40 

260 

11.00 

2.89 

-131.9 

197.0 

0.82 

130 X 130 

1.4 

g278 

277.9 

0.0 

40 

990 

40 

260 

11.00 

2.89 

-99.8 

229.2 

0.82 

130 X 130 

1.3 

g288 

287.9 

0.0 

40 

990 

40 

260 

11.00 

2.89 

-112.1 

216.8 

0.82 

130 X 130 

1.6 

g298 

297.9 

0.0 

40 

990 

40 

260 

11.00 

2.89 

-132.7 

196.2 

0.82 

130 X 130 

1.4 

g308 

307.9 

0.0 

40 

990 

40 

260 

11.00 

2.89 

-145.1 

183.9 

0.82 

130 X 130 

1.5 

g318 

317.9 

0.0 

40 

990 

40 

260 

11.00 

2.89 

-161.6 

167.4 

0.82 

130 X 130 

1.6 

g328 

327.9 

0.0 

40 

990 

40 

260 

11.00 

2.89 

-186.3 

142.6 

0.82 

130 X 130 

1.7 

g338 

337.9 

0.0 

40 

990 

40 

260 

11.00 

2.89 

-202.8 

126.1 

0.82 

130 X 130 

1.9 

g348 

347.9 

0.0 

40 

990 

40 

260 

11.00 

2.89 

-247.3 

164.9 

0.82 

130 X 130 

1.9 

g353 

353.0 

0.0 

40 

990 

40 

260 

11.00 

2.89 

-300.1 

150.0 

0.82 

145 X 90 

2.6 

gOlO 

10.0 

0.0 

40 

853 

40 

260 

9.48 

2.89 

-145.1 

266.3 

0.82 

160 X 110 

1.8 

g015 

16.8 

0.0 

40 

660 

40 

260 

7.33 

2.89 

-128.6 

200.3 

0.82 

200 X 130 

1.4 



Table 4. Properties of the ATCA and Parkes combined H l cubes containing continuum emission. Column headings 

are the same as in Table 3. 


cube 

1 

n 

b 

n 

61 

n 

npixi 

6b 

n 

npixb 

Al 

n 

<1 

'^min 

(km s“^) 

Vmax 

(km s“^) 

6v 

(km s“^) 

9 

{" X ") 

rms 

(K) 

g255 

255.4 

0.0 

-0 

721 

30 

351 

-6.0 

2.9 

-227.5 

265.5 

0.82 

143 X 110 

1.2 

g260 

260.4 

0.0 

30 

721 

30 

351 

-6.0 

2.9 

-227.5 

265.5 

0.82 

124 X 120 

1.4 

g265 

265.5 

0.0 

30 

721 

30 

351 

-6.0 

2.9 

-227.5 

265.5 

0.82 

138 X 110 

1.8 

g270 

270.4 

0.0 

30 

721 

30 

361 

-6.0 

3.0 

-227.5 

265.5 

0.82 

132 X 114 

1.1 

g275 

275.4 

0.0 

30 

721 

30 

361 

-6.0 

3.0 

-227.5 

265.5 

0.82 

131 X 116 

1.2 

g280 

280.4 

0.0 

30 

721 

30 

346 

6.0 

2.9 

-227.5 

265.5 

0.82 

133 X 112 

1.2 

g285 

285.4 

0.0 

30 

721 

30 

351 

6.0 

2.9 

-227.5 

265.5 

0.82 

126 X 115 

1.8 

g290 

290.4 

0.0 

30 

721 

30 

351 

6.0 

2.9 

-227.5 

265.5 

0.82 

129 X 112 

1.5 

g295 

295.4 

0.0 

30 

721 

30 

351 

6.0 

2.9 

-227.5 

265.5 

0.82 

123 X 116 

1.4 

gSOO 

300.4 

0.0 

30 

721 

30 

351 

6.0 

2.9 

-227.5 

265.5 

0.82 

124 X 113 

1.6 

g305 

305.4 

0.0 

30 

721 

30 

351 

6.0 

2.9 

-227.5 

265.5 

0.82 

135 X 112 

1.3 

g310 

310.4 

0.0 

30 

717 

30 

348 

6.0 

2.9 

-227.5 

265.5 

0.82 

131 X 119 

1.8 

g315 

315.4 

0.0 

30 

716 

30 

346 

6.0 

2.9 

-227.5 

265.5 

0.82 

131 X 116 

1.7 

g320 

320.4 

0.0 

30 

721 

30 

346 

6.0 

2.9 

-227.5 

265.5 

0.82 

129 X 119 

1.7 

g325 

325.4 

0.0 

30 

721 

30 

351 

6.0 

2.9 

-227.5 

265.5 

0.82 

142 X 117 

1.7 

g330 

330.4 

0.0 

30 

721 

30 

356 

6.0 

3.0 

-227.5 

265.5 

0.82 

129 X 115 

2.2 

g335 

335.4 

0.0 

30 

716 

30 

351 

6.0 

2.9 

-227.5 

265.5 

0.82 

137 X 126 

2.2 

g340 

340.4 

0.0 

30 

716 

30 

351 

6.0 

2.9 

-227.5 

265.5 

0.82 

127 X 121 

2.1 

g345 

345.4 

0.0 

30 

716 

30 

346 

6.0 

2.9 

-227.5 

265.5 

0.82 

147 X 114 

2.6 

g350 

350.4 

0.0 

30 

716 

30 

351 

6.0 

2.4 

-227.5 

265.5 

0.82 

152 X 108 

2.4 

g355 

355.6 

0.0 

30 

761 

30 

351 

6.3 

2.9 

-227.5 

265.5 

0.82 

161 X 108 

3.4 

gOlO 

9.5 

0.0 

30 

1200 

30 

350 

10.0 

2.9 

-227.5 

265.5 

0.82 

190 X 93 

2.0 

g015 

16.5 

0.0 

30 

860 

30 

350 

7.2 

2.9 

-227.5 

265.5 

0.82 

210 X 92 

1.9 
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I g295_045 1.4219 GHz 



-2 0 2 
u (kA) 


Fig. 1.— The u-v coverage for a typical SGPS I pointing in the held g295. The 6 km antenna 
has been omitted from this plot to show the coverage of the inner u-v plane. These data use 
only the east-west conhgurations of the ATCA. 


I g012_045 1.4220 GHz 



-1.5 -1 -0.5 0 0.5 1 1.5 

u (kX) 


Fig. 2.— The u-v coverage for a typical SGPS II pointing in the held g012. The 6 km 
antenna has been omitted from this plot to show the coverage of the inner u-v plane. The 
u-v coverage of this pointing dihers from the coverage shown in Figure 1 because of the use 
of the north-south hybrid arrays in SGPS II and the lower declination of the held. 
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Parkes Data 255 de 9 < I < 318 deg 
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Parkes Data 318 deg < I < 21 deg 


54.405958 km/s 
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Fig. 3.— Parkes data at u = 34.4 km s“^ for the region 255° < I < 318° (a) and 318° < 
I < 21° (b). The greyscale is linear from 3 to 100 K and has been intentionally saturated to 
bring out the low-level emission. 


SGPS II cubes 


SGPS I cubes 



Galactic Longitude (deg) Galactic Longitude (deg) 

Fig. 4.— Coverage of the SGPS I & II continuum subtracted combined cubes. Field names 
are the names used for the data cubes as released and as given in Table 3. 
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Galactic Longitude 


Fig. 5.— Longitude-velocity diagram of continuum subtracted data for the SGPS I region 
at 6 = 0°. The greyscale runs from 3 to 120 K and uses a square-root transfer function. 
The occasional white vertical stripes (e.g. near I = 306°) in the image are the result of 
H I absorption towards strong continuum sources that have been subtracted. Similarly, the 
horizontal stripes, especially near n = 0 km s“^ are due to H l self-absorption in the local 


gas. 
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Fig. 6.— Longitude-velocity diagram of continuum subtracted data for the SGPS II region 
at 6 = 0°. The greyscale runs from 3 to 120 K and uses a square-root transfer function. As 
in Figure 5, the vertical stripes seen here are due to continuum absorption and the horizontal 
stripes near n = 0 km s“^ are due to H l self-absorption. 
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Velocity: 75.S5 km/s 
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Fig. 7.— Comparison of Parkes and Combined data for the g268 cube at u = 75.05 km s“^. 
The beam size is marked as an open circle in the upper left-hand corner of each image. The 
greyscale is the same for the two panels, running linearly from 1 to 75 K. 



